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Abstract 

The J1-J2 model with the biquadratic (plaquette-four-spin) interaction was simulated with the 
numerical-diagonalization method. Some limiting cases of this model have been investigated thor- 
oughly. Taking the advantage of the extended parameter space, we survey the phase boundary 
separating the Neel and valence-bond-solid phases. According to the deconfined-criticality scenario, 
the singularity of this phase boundary is continuous, accompanied with unconventional critical in- 
dices. Diagonalizing the finite-size cluster with < 36 spins, we observe a signature of continuous 
phase transition. Our tentative estimate for the correlation- length critical exponent is = 1.1(3). 
In order to elucidate a non-local character of criticality, we evaluated the Roomany-Wyld f3 function 
around the critical point. 
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I. INTRODUCTION 



The deconfined criticalityi"- is arousing much attention recently.- Naively,- the (ground- 
state) phase transition separating the Neel and valence-bond-solid (VBS) phases is dis- 
continuous, because the adjacent phases possess distinctive order parameters such as the 
sublattice magnetization and the dimer- coverage pattern, respectively. However, according 
to the deconfined-criticality scenario, the transition is continuous, accompanied with novel 
critical indices;-""— afterward, we make an overview on recent computer-simulation studies. 

In this paper, we investigate an extended version of the J1-J2 model^"— [see Eq. ([T])] 
by means of the numerical diagonalization method for the cluster with < 36 spins. The 
extension of the parameter space, {J2,Q), permits us to investigate the above-mentioned 
criticality from a global viewpoint. A schematic phase diagram is presented in Fig. [H 
details are explicated afterward. As mentioned above, the critical branch (dashed line) is 
our concern. 

The Hamiltonian of the J1-J2 model with the biquadratic interaction Q is given by 
n = Ji J]Si-S,+J2 J2 S.-S,-Q J][(S,-S,-l/4)(Sfc-S^-l/4)+(SrS^-l/4)(S,-Sfc-l/4)]. 

(ij) m) mi] 

(1) 

Here, the quantum spin-1/2 operators {S,} are placed at each square-lattice point i. The 
summations, ^[ijki]^ '^^^^ possible nearest-neighbor pairs {ij), next- 

nearest-neighbor pairs {{ij)), and plaquette spins [ijkl], respectively; here, the arrangement 
of indices, [ijkl], around a plaquette, □, is ^D^. The parameters, Ji, J2, and Q, are the 
corresponding coupling constants. Hereafter, we consider Ji as the unit of energy; namely, 
we set Ji = 1. Both frustration and biquadratic interactions, J2 and Q, stabilize the VBS 
phase. 

In some limiting cases (subspaces), namely, either Q = or J2 = 0, the Hamiltonian ([1]) 
has been investigated extensively so far. The case of Q = 0, namely, the J1-J2 model,— 
has been analyzed with the series-expansion methocii^""— and the numerical diagonaliza- 
tion method.—"— (The quantum Monte Carlo method is inapplicable to such a frustrated 
magnetism because of the negative-sign problem.) The Neel, VBS, and coUinear phases 
appear— successively in the regimes, J2/J1 ^ 0.4, 0.4 < J2/J1 < 0.6, and 0.6 < J2/ Ji, 
respectively. The singularity-^i at J2/J1 ~ 0.4 (0.6) is continuous (discontinuous^^). On the 
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one hand, the case of J2 = 0, namely, the J-Q model^^i^^ {Ji-Q model in our notation), 
admits the use of the quantum Monte Carlo method. The Neel phase turns into the VBS 
phase at a considerably large biquadratic interaction Q/J ^ 25.^^ This fact indicates that 
the antiferromagnetic order at J2 = is considerably robust against Q. Nevertheless, the 
large-scale quantum Monte Carlo simulation admits detailed analyses of criticality. The 
correlation-length critical exponent is estimated as z/ = 0.78(3) (Ref. [2^ and u = 0.68(4) 
(Ref. I25I ). According to Ref. |26!, logarithmic scaling corrections prevent us from taking 
the thermodynamic limit reliably. Possibly because of this difficulty, it is still controver- 
sial whether the singularity is continuous or belongs to a discontinuous one accompanied 
with an appreciable latent heat^~— ; see aforementioned Ref. as well. Last, it has be 
mentioned that there have been reported a number of intriguing extentions such as the 
third-neighbor interaction,— internal symmetry enlargement,— and spatial anisotropy.— "— 
The dimer model^i^ is also a clue to the study of deconfined criticality. 

In this paper, taking the advantage of the extended parameter space (J2,<5) (see Fig. 
[1]), we survey the critical branch in details. Scanning the parameter space, we found that 
around an intermediate regime Q = 0.2, the finite-size-scaling behavior improves signifi- 
cantly. Around this optimal regime, we carry out finite-size-scaling analyses 

The rest of this paper is organized as follows. In Sec. [Tll we present the simulation results. 
Technical details are explicated. In Sec. IIIIl we address the summary and discussions. 



II. NUMERICAL RESULTS 

In this section, we present the numerical results for the J1-J2-Q model, Eq. ([T]). We 
employ the exact-diagonalization method for the cluster with A^(= L^) < 6^ spins. Our aim 
is to clarify the nature of the phase boundary separating the Neel and VBS phases. For that 
purpose, we scrutinize the excitation gap 

AE = -E'[(o,7r),o,+,+,+] ~ -E'[(o,o),o,+,+,+]- (2) 

Here, the variable Ea denotes the ground-state energy within the sector (subspace) spec- 
ified by the set of quantum numbers, a = {k, S^,±,±,±). These quantum numbers 
{k, S^, ±, ±, ±) are associated with the symmetry groups such as translation along the 
rectangular edges, internal-spin 2;-axis rotation (total-spin conservation along the z axis), 
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internal-spin inversion (5*/ f-)- —S^), 9 = n lattice rotation, and lattice inversion, respec- 
tively. These symmetry groups commute mutually, because we restrict ourselves within the 
sector with k^^ = and = 0. The excitation gap, Eq. ([2]), opens (closes) for the Neel 
(VBS) phase in the thermodynamic limit. (The degeneracy in the VBS phase corresponds to 
the translationally- invariant dimer-coverage patterns.) The excitation branch ([2]) does not 
correspond to the first energy gap: There appear numerous low-lying excitations due to the 
spontaneous-symmetry breakings in both Neel and VBS phases. Here, aiming to discrimi- 
nate these phases clearly, we look into an excitation branch ([2]), which is a key ingredient of 
the finite-size-scaling analysis. 

A. Survey of the phase boundary between the VBS and Neel phases 

In Fig. m we plot the scaled energy gap LAE for Q = 0.2, various J2, and = 4, 
16, and 36. Here, the variable L{= \/N) denotes the linear dimension of the cluster. The 
intersection point of the curves indicates a location of the critical point. We see that a 
transition occurs at J2 ~ 0.3. 

In order to estimate the transition point precisely, in Fig. [3l we plot the approximate 
critical point J2c{Li, L2) for 2/(Li + L2) with 2 < Li < L2 < 6; the parameters are the same 
as those of Fig. [21 Here, the approximate critical point J2c{Li, L2) denotes a scale-invariant 
point with respect to a pair of system sizes (Li,L2). Namely, it satisfies the relation 

-^l^-^(-^l)|j2 = J2c(Ll,L2) = L2AE{L2)\j^ = J2,{Li,L2)- (3) 

The least-squares fit to the data of Fig. [3] yields an estimate J2C = 0.273(12) in the ther- 
modynamic limit L —7- 00. Replacing the abscissa scale with l/L"^, we obtain an alternative 
estimate J2C = 0.3095(91). The discrepancy between these estimates, ~ 3.7-10"^, appears to 
dominate the insystematic (statistical) error, ^ 10^^. Considering the former as an indicator 
of the error margin, we estimate the transition point J2C = 0.273(37). 

We carried out similar analyses for various values of Q. The results are shown in Fig. 
m The error margin is suppressed around an intermediate regime Q ~ 0.1. As a matter of 
fact, in Fig. |5l we show that the scaling behavior becomes optimal around Q ~ 0.2. (Such 
a redundancy is a benefit of the parameter-space extension.) On the contrary, the error 
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margins get enhanced around the multi-critical point as well as in the large-Q regime. 
As mentioned in the Introduction, the limiting cases, Q = and J2 = 0, have been 



studied extensively.— 1^ In Ref. 
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The transition point ( J2, Q) ~ (0.4, 0) was reported. 
This result agrees with ours. On the contrary, the transition point^*^ ( J2, Q) ~ (0, 25) is not 
reproduced by the present approach. This discrepancy may be attributed to the logarithmic 
corrections^^ inherent in the subspace J2 = 0; these corrections prevent us from taking the 
thermodynamic limit reliably. Around Q ~ 0, there arise pronounced scaling corrections, 
which cannot be appreciated properly; this regime would not be accessible by the numerical 
diagonalization method. As a matter of fact, the Q-driven VBS state is so unstable that 
a considerably weak antiferromagnetic interaction Ji = 0.04 destroys VBS immediately. 
This difficulty is remedied by the inclusion of the magnetic frustration J2, which appears 
to stabilize VBS significantly (Fig. H]), and recover an applicability of the diagonalization 
method. 

Last, we make a comment. The present simulation cannot rule out a possibility of the 
discontinuous phase transition. Even in such a case, the phase boundary determined through 
the criterion, Eq. ([3]), makes sense, albeit the value of LAE at the transition point no longer 
converges to a universal value (critical amplitude). 



B. Correlation-length critical exponent 

In this section, we estimate the correlation-length critical exponent u. In Fig. 0, we plot 
the approximate critical exponent 

,^ . ^ ln(Li/L2) , . 

^' ^^-ln{9,[LiAE(Li)]/9,[L2AE(L2)]}U,=.,,(z.,,,)' 

for 2/{Li + L2) and 2 < Li < L2 < 6. The parameters are the same as those of Fig. [2l The 
least-squares fit to these data yields u = 1.094(26) in the thermodynamic limit. In order to 
appreciate possible systematic errors, we made an alternative extrapolation with the 1 / L^- 
abscissa scale. Thereby, we arrive at u = 1.087(14). In this (optimal) parameter space, 
Q = 0.2, the systematic error (deviation between different abscissa scales) appears to be 
negligible. On the one hand, at Q = 0.3, we estimated z/ = 0.901(21) and 0.988(17) through 
the 1/L- and 1/L^-extrapolation schemes, respectively. The discrepancy (systematic error) 
appears to be enhanced by the scaling (possibly, logarithmic-2&) corrections. At Q = 0.1, 



via the 1/L- and 1/L^-extrapolation schemes, we arrive at z/ = 1.316(87) and 1.198(55), 
respectively. Both systematic and insystematic errors get enhanced in the proximity to the 
multicritical point. Considering these deviations of various characters as an indicator of the 
error margin, we estimate the critical exponent as 

z/ = 1.1(3). (5) 

A validity of the scaling analysis is considered in the next section. 



C. Non-local behavior of the /3 function 

In Fig. El we plot the (3 function, l3i'^^{J2)^ for Q = 0.2, various J2, and (^1,^2) = (+) 
(2, 4), (x) (4, 6), and (*) (2, 6). Here, we calculated the /3 function with the Roomany-Wyld 
f o rixxu. 1 3)"" 

nRw rj ^ ^ l + HAE{L,)/AE{L,))/\n{L,/L2) 

for a pair of system sizes (Li, L2). 

The /3 function elucidates a non-local character of the criticality. In the vicinity of the 
critical point J2C, the /3 function falls into a linearized formula 

/3(J2) = ^(J2c- J2), (7) 

with the correlation-length critical exponent v. In Fig. |6l we present the theoretical pre- 
diction (3 = (0.273 — J)/l.l (dotted line); here, the critical parameters, J2C = 0.273 and 



1.1, are estimated in the preceding sections. The behavior of (3]^^^ is consistent with 
the theoretical prediction for a considerably wide range of J2, validating out treatment. 
In particular, the linearity of Pl^^ is remarkable, suggesting that scaling corrections are 
eliminated satisfactorily, at least, around Q = 0.2. Possibly, the J2-driven phase transition, 
rather than the Q-driven transition, is reasonable in the sense of the renormalization-group 
flow. 



D. Ground-state energy: An indication of the VBS-collinear phase transition 

In Fig. [TJ we present the ground-state energy per unit cell, Eg/N, with = 36 for Q = 0.2 
and various J2- From the figure, we observe a cusp-like anomaly around J2 ~ 0.7. This 
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anomaly indicates an onset of the first-order phase transition between the VBS and colhnear 
phases. Making similar analyses for various values of Q, we obtained a schematic feature of 
the VBS-coUinear phase boundary as shown in Fig. [T] Character of the multicritical point 
is unclear. 



III. 



SUMMARY AND DISCUSSIONS 



The phase diagram of the extended J1-J2 model with the biquadratic interaction Q, Eq. 
(II]), was investigated by means of the numerical diagonalization method for N < 36. Taking 
the advantage of the extended parameter space, {J2,Q), we surveyed the critical branch 
separating the Neel and VBS phases (Fig. H]) in particular. The finite-size corrections 
appear to be suppressed around an intermediate regime Q ~ 0.2. Encouraged by this 
finding, we analyzed the criticality around Q ^ 0.2 with the finite-size-scaling method. As a 
result, we estimate the correlation- length critical exponent z/ = 1.1(3). A non-local feature 
of the (3 function is accordant with the theoretical prediction (Fig. |6]), giving support to our 
treatment. 

As a reference, we overview related studies. For the J-Q model (square-lattice anti- 



ferromagnet with t 
u = 0.68(4) (Ref. 



le biquadratic interaction), the estimates, u = 0.78(3) (Ref 



and 



25l ). were rep orted. As for the spatially-anisotropic spin systems, the 



exponent z/ = 0.80(15) (Ref. 



Sa) and z/ = 0.92 



10) (Ref. |39i) were obtained. These results 



are to be compared with u = 0.7112(5) (Ref. l43l ) for the d = 3 Heisenberg universality class. 
Our result indicates a tendency toward an enhancement for the correlation-length critical 
exponent, as compared to that for the d = 3 Heisenberg universality class. 

According to Ref. information on the deconfined criticality can be drawn from the 
details of the multi-critical point. This viewpoint would provide an evidence whether the 
criticality belongs to the deconfined criticality; at present, it is not conclusive whether the 
Neel- VBS transition belongs to the deconfined criticality. This problem will be addressed in 
future study. 
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FIG. 1. A schematic (ground-state) phase diagram for the J1-J2-Q model ([T]) is presented. The 
sohd (dashed) hne stands for the phase boundary of discontinuous (continuous) character. We 
investigate the phase boundary separating the Neel and VBS phases. Details of the multi-critical 
point and the VBS-collinear phase boundary are uncertain. 




FIG. 2. The scaled energy gap LAE is plotted for the next-nearest-neighbor interaction J2 and 
the system size L =(+) 2, (x) 4, and (*) 6. The biquadratic interaction is fixed to Q = 0.2. The 
data indicate an onset of continuous phase transition around J2 ~ 0.3. 
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FIG. 3. The approximate critical point J2ciLi,L2) ([3]) is plotted for 2/(Li + L2) with 2 < Li < 
L2 < 6. The parameters are the same as those of Fig. El The least-squares fit to these data yields 
J2c = 0.273(12) in the thermodynamic limit. A possible systematic error is considered in the text. 
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FIG. 4. The phase boundary between the Neel and VBS phases (critical branch) is presented. 

4 R. R. P. Singh, Physics 3, 35 (2010). 

5 O. I. Motrunich and A. Vishwanath, Phys. Rev. B 70, 075104 (2004). 

6 A. Tanaka and X. Hu, Phys. Rev. B 74, 140407(R) (2006). 
G.-Z. Liu, Phys. Rev. B 71, 172501 (2005). 

^ R. Dillenschneider and J. Richert, Phys. Rev. B 73, 224443 (2006). 

9 P. Ghaemi and T. Senthil, Phys. Rev. B 73, 054415 (2006). 

° I. O. Thomas and S. Hands, Phys. Rev. B 75, 134516 (2007). 

1 D. H. Kim, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. 79, 2109 (1997). 

2 M.P. Gelfand, R.R.P. Singh, and D.A. Huse, Phys. Rev. B 40, 10801 (1989). 
^ E. Dagotto and A. Moreo, Phys. Rev. Lett. 63, 2148 (1989). 

^ H. J. Schulz and T. A. L. Ziman, Europhy. Lett. 18, 355 (1992). 

^ J. Oitmaa and Z. Weihong, Phys. Rev. B 54, 3022 (1996). 

9 



1.12 





1.11 - 


CM 


1.1 - 




1.09 - 


-J 


> 


1.08 - 




1.07 - 




1.06 - 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 

2/(Li+L2) 



FIG. 5. The approximate critical exponent v{Li^L2) dH) is plotted for 2/(Li + L2) with 2 < Li < 
L2 < 6. The parameters are the same as those of Fig. El The least-squares fit to these data yields 
V = 1.094(26) in the thermodynamic limit. A possible systematic error is considered in the text. 
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a wide range of J2. 
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